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ABSTRACT 


It  is  considered  likely  that  a  meteorological  rocket  which  would  have 
the  brief  capability  of  a  meteorological  satelite  will  soon  be  needed.  Such 
a  vehicle  could  be  used  to  test  out  new  ideas  and  equipment  for  meteoro¬ 
logical  satellites  and  to  provide  a  quick  satellite's  eye  view  of  a  question¬ 
able  weather  area. 

It  will  be  necessary  for  this  type  of  vehicle  to  contain  some  type  of 
attitude  control.  This  report  describes  the  factors  involved  in  attitude 
control  and  develops  a  theoretical  solution.  Various  known  values  are 
then  used  with  the  formulas  proposed  to  demonstrate  the  feasibility  of 
the  solution.  It  is  concluded  that  the  problem  could  be  solved  by  applying 
implusive  torquing  to  bring  the  longitudinal  axis  of  the  rocket  into  correct 
alignment  to  produce  the  desired  attitude.  The  torque  could  be  applied 
by  a  cold  propellant  discharged  by  a  bang-bang  servo  mechanism. 
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SECTION  1  -  INTRODUCTION 


Meteorological  satellites,  such  as  Tiros,  are  proving  invaluable  as 
an  eye  in  the  sky  to  observe  general  weather  patterns;  however,  small 
vertical  probe  rockets,  such  as  L>oki  and  Areas,  are  necessary  for 
sampling  the  atmosphere  up  to  250,  000  ft.  It  is  considered  likely  that 
an  intermediate  rocket  (suborbital  vehicle)  will  soon  be  necessary.  Such 
a  vehicle  would  have  the  brief  capability  of  a  meteorological  satellite, 
but  the  cost  and  mobility  of  a  vertical  probe. 

Two  possible  uses  of  such  a  vehicle  would  be: 

1.  It  could  be  used  to  test  out  new  ideas  and  equipment  for 
meteorological  satellites. 

2.  It  could  be  used  to  give  a  quick,  satellite's  eye  view  of  a 
questionable  weather  area,  such  as  a  region  where  there  is  an  alleged 
tropical  cyclogenesis.  An  existing  satellite  might  have  to  make  many 
orbits  before  it  would  be  in  a  position  that  a  meteorological  rocket  could 
assume  almost  immediately. 

Most  of  the  design  of  an  intermediate  vehicle  will  probably  be  modi¬ 
fications  of  larger  or  smaller  types  of  rockets;  however,  it  will  be  nec¬ 
essary  for  the  intermediate  vehicle  to  have  some  means  of  attitude  con¬ 
trol.  The  problem  of  attitude  control  for  these  vehicles  will  be  unique. 

The  purpose  of  this  report  is  to  present  a  theoretical  solution  which 
can  be  used  as  a  basis  for  the  development  of  a  low-cost  attitude  control 
system.  The  solution  to  the  problem  is  approached  in  the  following 
phases: 


1.  Describe  the  attitude  control  problem. 

2.  Define  a  mathematical  system  and  symbology  that  can  be 
used  to  solve  the  problem. 

3.  State  the  equations  of  angular  motion  for  a  body  in  free  fall. 

4.  Develop  a  solution  for  the  problem. 
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5.  Examine  some  of  the  state-of-the-art  values  to  determine 
if  the  solution  approach  is  practical. 

The  solution  to  the  problem  of  attitude  control  presented  in  this  re¬ 
port  was  developed  by  the  author,  based  on  ideas  on  the  general  subject 
initiated  by  Capt  J.  H.  Jacobsmeyer,  USAF,  and  Dr.  R.  M.  Howe, 
Professor  of  Aeronautics  and  Astronautics,  University  of  Michigan. 


SECTION  2  -  DESCRIPTION  OF  THE  PROBLEM 


A  vertical  or  horizont.^l  attitude  is  the  most  likely  attitude  that  the 
payload  would  be  required  to  assume.  A  vertical  attitude  is  defined  as 
an  attitude  in  which  the  longitudinal  axis  of  the  payload  is  normal  to  the 
plane  of  the  horizon.  A  horizontal  attitude  is  defined  as  an  attitude  in 
which  the  longitudinal  axis  of  the  payload  is  in  a  horizontal  plane  and 
also  in  the  plane  of  the  trajectory. 

In  addition,  it  will  probably  be  necessary  to  spin  the  vehicle  about 
its  longitudinal  axis  for  reasons  of  stability.  The  sensor  used  in  this 
vehicle  will  probably  be  some  type  of  detector  such  as  a  radiometer  or 
a  camera. 

Fig.  1  shows  a  vertically-oriented  payload.  If  the  detector  is  aimed 
out  of  the  side  of  the  payload  at  a  constant  angle  ({>,  and  if  the  trajectory 
does  not  carry  it  significantly  downrange,  the  detector  will  scan  the  earth 
below  in  nearly  a  circle.  If  angle  (J>  is  allowed  to  vary  uniformly,  the 
scanner  will  describe  a  spiral.  Since  the  vehicle's  actual  trajectory  will 
be  a  parabola  with  velocity  components  in  the  horizontal  and  vertical,  the 
resultant  detector  sweep  will  be  at  best  a  modified  prolate  cycloid. 

Fig.  2  shows  the  horizontally-oriented  payload.  In  this  case,  the 
detector  will  be  aimed  normal  to  the  spin  axis.  Thus  the  sweep  of  the 
detector  will  be  a  helix  as  the  vehicle  moves  downrange.  For  over  half 
of  each  revolution,  the  detector  will  be  directed  toward  space  leaving  the 
useful  portion  of  each  sweep  to  be  the  horizon-to-horizon  scan  of  the 
earth.  Each  successive  scan  of  the  earth  will  be  a  strip  farther  down- 
range.  The  overlap  or  gap  between  successive  scans  will  be  determined 
by  the  height,  downrange  velocity  component,  detector  beam  width,  and 
spin  rate. 

The  attitude  control  problem  of  reaching  the  horizontal  is  basically 
the  same  as  that  of  reaching  the  vertical;  however,  the  equipment  nec¬ 
essary  to  reach  and  hold  the  horizontal  varies  from  that  required  to 
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Fig.  1:  Vertically-Oriented  Payload. 


reach  and  hold  the  vertical.  In  reaching  the  vertical,  the  horizontal 
plane  can  be  easily  sensed  by  a  spinning  radiometer;  therefore,  if  the 
horizontal  plane  is  known,  the  payload's  vertical  is  uniquely  defined.  If 
a  payload  were  required  to  assume  an  attitude  other  than  the  vertical,  a 
gyro  or  celestial  navigation  system  would  be  required.  A  difference, 
then,  exists  in  cost,  as  a  radiometer  is  inexpensive  as  compared  to  a 
gyro  or  celestial  navigation  system. 


3 


APGC-TDR-62-30 


Fig.  2:  Horizontally-Oriented  Payload. 


Since  bringing  the  payload  to  the  horizontal  is  the  more  general 
case,  it  will  be  discussed  exclusively  hereafter. 

SECTION  3  -  DEFINITION  OF  TERMS 

BODY -FIXED  COORDINATE  SYSTEM 

Most  components  will  be  taken  in  a  cartesian,  body-fixed,  principal 
axis  system  using  axes  which  are  called  x,  y,  and  z  (see  Fig.  3).  The 
origin  is  at  the  center  of  mass. 

z  -  the  axis  of  mass  symmetry  of  the  payload.  It  is  the  longitu¬ 
dinal  axis. 
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Fig.  3:  Body-Fixed  Coordinate  System. 


X  and  y  -  the  lateral  .axes  passing  through  the  center  of  mass  and 
rotating  with  the  body. 


TRAJECTORY  COORDINATES 

These  axes  are  determined  by  the  earth  and  by  the  trajectory.  They 
also  are  cartesian  but  nonrotating.  Their  origin  is  also  the  center  of 
mass,  with  the  axes  named  u,  v,  and  w  (see  Fig  4). 

u  -  this  axis  is  horizontal  and  in  the  plane  of  the  trajectory.  The 
term  horizontal  axis  will  be  used  in  this  report  to  indicate  the  u  axis. 

V  -  the  vertical  axis  is  defined  as  being  normal  to  the  horizontal 
plane. 

w  -  this  axis  is  horizontal  but  perpendicular  to  the  plane  of  the 
trajectory.  The  term  normal  axis  will  be  used  in  this  report  to  indicate 
the  w  axis. 


5 


APGC-TDR-62-30 


Fig.  4:  Trajectory  Coordinates. 


SUBSCRIPTS 

Subscripts  will  designate  along  which  axis  the  components  are  taken. 


DIFFERENTIAL  OPERATOR 

The  term  p,  as  used  in  this  report  is  the  same  as  the  time  derivative 
often  seen  as  (d/dt)  or  (*).  The  use  of  p  is  for  simplicity  and  convenience 
in  converting  into  Laplace  transforms. 

For  example,  pM  is  identically  the  same  as  (dM/dt)  or  (M)  or  the 

first  time  derivative  of  M,  while  p^M  is  the  same  as  (d^M/dt^)  or  (M) 
or  the  second  time  derivative  of  M. 
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VECTORS 

M  -  the  total  external  moment  or  torque  acting  on  a  body. 

H  -  the  total  angular  momentum  of  a  body, 
n  -  the  total  angular  velocity  of  a  body. 

R  -  the  position  vector  from  the  center  of  mass  to  the  z  axis  linear 
accelerometer. 

r  -  the  position  vector  from  the  center  of  mass  to  the  control  jet. 
a  -  linear  acceleration. 


ANGULAR  MEASUREMENTS 

6  -  the  angle  between  the  u  axis  and  the  z  axis. 

-p6  -  the  closing  r.ate  of  the  u  axi  3  and  the  z  axis. 


TENSOR 

I  -  the  inertia  tensor. 


COMPONENTS 


Since  we  chose  principal  axes,  the  products  of  inertia  in  the  inertia 
tensor  become  zero. 

^x  "^x  "^zx 

1  =  J  1  J 

xy  y  zy 

•^z 

Since  vectors  M,  H,  and  fl  are  three  dimensional,  they  are  written 
in  matrix  form  as  follows 


X 

O 

O 


O 

I 


O 

o 


y 

o  i 


diagonal  matrix. 


m 

h 

(a) 

X 

X 

X 

,  H  = 

h 

,  a  = 

CO 

y 
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y 
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co„  • 

z 
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is  the  component  of  Q  in  the  x-y  plane,  i.  e. 
c-y  ^~y  *  y 


SCALARS 

t  -  time. 

k  -  linear  torque  coefficient, 
f  -  the  viscous  damping  coefficient. 


SECTION  4  -  mathematical  DESCRIPTION 

.  .  *  *  '  t 

In  any  body,  the  sum  of  the  torques  equals  the  rate  of  change  of  the 
angular  momentum.  The  equation  is  . 

M  =  pH.  (1) 

Vector  equation  1  may  be  broken  up  into  components.  If  derivatives 
are  taken  in  a  rotating  axis  system,  there  must  be  coriolis  terms,  so 
that  equation  1  becomes 


mx  =  phx 

-  hyWjJ 

+  ^7.^} 

m  =  ph 

-  h  w 

+  h  w 

y  y 

Z  X 

X  2 

m  =  ph 

-  h  w 

+  h  w 

z  ^  z 

X  y 

y 

The  definition  of  angular  momentum  is 

H  =  in 


where  H  is  a  column  matrix  (h^^,  hy,  h^);  I  is  a  diagonal  matrix 
(Ijj,  ly,  I^):  and  n  is  a  column  matrix  Wy,  w^).  Since  x,  y,  and  z 

are  principal  axes,  the  components  produce  three  separate  equations: 

hj^  —  (4) 

hy  =  lyWy 

^Z  ”  ^z***  z , 
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There  is  mass  symmetry  about  the  z  axis,  so  we  can  say 

1=1.  (5) 

X  y  '  ' 

If  we  consider  a  case  where  no  torques  are  applied  about  the  z  axis, 
equation  2  reduces  to 

"'x  =  +  <Iz  -  ^>“z“y  (6a) 

"V  =  ^xP“y  ■  <^z  ■  Ix>  “z“x-  (6b) 

If  we  assume  that  the  rocket  is  spinning  about  the  z  axis  at  a  con¬ 
stant  rate  and  if  the  applied  torques  are  known,  then  equations  6a 
and  6b  are  coupled,  first  order,  differential  equations.  These  coupled 
equations  combine  to  give  for  the  characteristic  equation: 

p^Wy  +  (  1  -  Iz/Ix)^“z ‘*»y  =  0*  (7b) 

In  the  case  of  the  meteorological  rocket,  Ij^  is  much  greater  than  I^, 
so  that  equations  7a  and  7b  reduce  to  approximately 

‘‘^x  =  ‘^z^  “x  "  ® 

p^Wy  +  w^Wy  =  (p^  +  w|)  Wy  r  0.  (7d) 

The  solutions  to  these  equations  are 


‘‘^x  "  ‘‘k-y 


sin  Ujjt 


‘**y  =  ‘*'x-y “z'^* 

The  meaning  of  these  equations  might  be  reasoned  physically.  That 
is  to  say,  most  of  the  total  angular  velocity  fl  is  contained  in  component 
w^.  The  remaining  angular  velocity  lies  in  the  x-y  plane,  or 

2  2  2  2  2  2 
0  -  ‘‘^x-y  "  “z  ‘^x  ‘^y.  (9) 


Since  axes  x  and  y  are  rotating  with  the  body,  components  and  co^ 

are  oscillating  sinusoidally  with  amplitudes  w  as  shown  by  equations 

X-y 

8a  and  8b.  Thus,  twice  during  each  revolution  Wj^  or  cOy  becomes  zero. 
At  the  instant  or  oj  becomes  zero,  the  moment  equations  6a' and  6b 
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become  uncoupled.  If  a  torque  (m  )  is  applied  in  an  impulsive  fashion 
each  time  is  zero,  the  payload  can  be  tilted  to  the  horizontal. 


SECTION  5  -  DESCRIPTION  OF  A  SIMPLE  MODEL 


The  attitude  control  problem  is  to  make  the  z  axis  of  the  payload 
coincide  with  the  u  axis  of  the  trajectory  as  quickly  as  possible,  and  to 
make  the  payload  maintain  this  attitude  during  the  flight. 

A  simple  two-degrees-of-freedom  gyro  can  be  used  to  maintain  the 
u  axis. 

A  bang-bang  control  jet  (cold  propellant)  can  be  used  to  bring  the  z 
axis  into  line.  The  control  jet  could  be  discharged  in  an  impulsive  man¬ 
ner  parallel  to  the  X  axis  at  distance  r  from  the  origin  (see  Fig.  2).  The 
jet  would  be  discharged  once  during  each  revolution  when  is  zero. 

The  integrated  effect  of  all  the  impulsive  torques  are  equivalent  to  a 
COatinuous  torque  (T). 

The  moment  of  inertia  ratio  (Ijj,/Ij^)  must  be  small  (0.05  or  less).  The 
spin  rate  must  also  be  small  (on  the  order  of  60  rpm).  If  the  torquing 
impulse  is  spread  over  1  radian  during  each  revolution,  the  imparted 
angular  momentum  (hy)  may  be  large  compared  to  the  original  angular 
momentum  (hj,).  If  the  impulses  are  applied  strong  enough  and  often 
enough,  the  torque  may  be  considered  continuous,  so  that  Newton's  Law 
of  Angular  Motion  may  apply  without  regard  to  small  nutation  and  pre¬ 
cession  effects:  therefore: 


T  =  I^p^  6.  (10) 

(Henceforth  I  will  be  referred  to  as  I.  ) 

It  is  apparent  from  equation  10  that  if  T  is  always  applied  in  a  direc¬ 
tion  to  drive  the  z  axis  toward  the  u  axis,  the  z  axis  will  oscillate  about 
the  u  axis  but  never  settle  on  it.  So  the  system  must  be  damped  by  a 
torque  (D)  opposite  to  p6. 

For  a  simplified  case,  the  system  would  be  linearized  and  ci  tically 
damped.  To  do  this,  T  would  be  made  proportional  to  the  6,  and  D  made 
proportional  to  p6. 
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6  can  be  measured  by  a  reference  to  the  attitude  gyro.  p6  can  be 
measured  by  a  linear  z-axis  acceleration,  since  the  only  significant  z- 
axis  acceleration  in  free  fall  will  result  from  the  centrifugal  force  pro¬ 
duced  by  the  angular  velocity. 

a^  =  j  p6  X  (Rxp6)  I  =  R  {p6)^.  (11) 

That  is,  if  a^  is  measured  by  the  accelerometer,  then  p6  is  also  deter¬ 
mined. 


To  make  T  proportional  to  6.  and  D  proportional  to  p6,  the  pulse 
length  or  pulse  strength  must  be  modulated.  This  can  be  done  simply 
only  over  a  rather  limited  dynamic  range. 


The  equation  of  motion  with  damping  thus  becomes 


m  =  T  +  D  =  I  p^  6. 

y 


(12) 


It  was  specified  that  m^  is  composed  of  T  and  D  whose  equivalents  are 


T  =  k6 


(13) 


-  D  =  fp6. 

So  that  equation  12  becomes  the  familiar  differential  equation 

Ip^6  +  fp6  +  k6  =  0. 


(14) 


(15) 


Normalizing  equation  I  5  gives 

p^6  +  (f/I)p6  +  (k/I)6  =  0.  (16) 

A  critically  damped  condition  is  one  in  which  there  is  no  overshoot 
nor  overdamping.  It  is  determined  by  the  condition  that 

f  =  4Ik.  (17) 


Thus  equation  15  becomes 

p^6  +  2  Vk/I  p6  +  (k/I)6  =  0.  (18) 

In  practical  engineering,  it  is  easier  to  set  up  equations  in  terms  of  the 
undamped  natural  frequency  (Wj^). 
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e 

where 


p  6  +  p6  +  6  =  0 


2 

cj  =  k/ 1, 
n 


(19) 


(20) 


To  see  if  the  above  equations  could  be  used  to  represent  a  workable 
system,  it  may  be  assumed  that  the  vehicle  would  have  a  moment  of 

inertia  <I)  of  about  1  60  slug  ft  .  This  moment  of  inertia  is  about  that  of 
the  J.avelin  vehicle  which  could  be  used  as  an  intermediate  meteorolog¬ 
ical  rocket. 


rocket  leaves  the  atmosphere,  its  inclination  angle  (6q) 
will  b  1.2  radians. 

For  a  first  approximation,  consider  a  critically  damped  case  where 
it  is  desired  to  tilt  from  an  initial  angle  of  70  degrees  (i.  e.  ,  6q  =  1-2 
radians)  to  within  5  degrees  of  the  horizontal  in  35  sec.  A  standard 
engineering  graph  gives  a  value  of  of  1/8  radian  per  second. 

From  equation  20.  k  may  be  found: 

2  222-2 
k  =  =  (160)(1/H)  slug  ft  rad  sec  =  2.  5  ft  lb.  (21) 

From  equation  17,  f  may  be  found  to  be  40  ft  lb  sec,  i.  e.  , 

f  =  ■|/4Ik  =  V4  X  160  X  2.  5  =  40  ft  lb  sec.  (22) 


Fig.  5:  History  of  6, 
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The  maximum  torque  which  might  be  expected  is 


Tmax  =  •'‘max  =  ^  ft  lb  .  1 . 2  rad  =  3  ft  lb.  (23) 

The  maximum  damping  torque  which  might  be  expected  is  about 

D  =  f(p6\  =  40  ft  lb  sec  x  0.  1  rad/sec  =  4  ft  lb.  (24) 

m^x  iTi3x 

For  the  first  approximation,  the  system  was  assumed  to  be  charac¬ 
terized  by  a  critically  damped  condition  from  the  beginning. 

As  a  second  approximation,  consider  a  more  practical  system  that 
is  initially  undamped.  When  6  reaches  a  predetermined  value 
=  0.7  radian)  damping  begins  (see  Fig.  5). 


The  history  of  the  torque  (T)  is  shown  as  a  function  of  6  in  Fig.  6. 
Fig.  7  is  a  graph  of  the  damping  torque  D.  Fig.  8  contains  complete 
graphical  representation. 


1.75 

1.50 

S  1.25 

-  1.00 

<1 

3 

J  0.75 
1  0.50 

o 

U 

0.25 

n 

— 

_ i 

B 

T,- 

1.75  ft/ 

lb 

m 

n 

_ _ J 

■ 

B 

y 

m 

■ 

Ql _ I _ I _ I _ I _ I - 1 - 1 - 1 

0  0.2  0.4  0.6  0.8  1.0  1.2  1.4  1.6 


S  (radians) 

Fig.  6:  Control  Torque  vs’ 6. 


APGC-TDR-62-30 


The  mathematics  of  the  system  are  as  follows: 

In  the  increment  0  to  1  (see  Fig.  8),  T  is  at  saturation  =  k6,  = 
2.  5  X  0.  7  =  1. 75  ft 


Tq  j  =  k6j  =  I(p  6)^  ^  =  saturation  torque, 


(25) 


so  that  the  angular  acceleration  in  the  region  is 

2  2  2 

(P  ^)q  j  “  /  ^  ft  lb/ 160  slug  ft  =  0,  01 1  rad/ sec  .  (26) 


With  this  constant  acceleration  it  takes  9.  5  sec  to  go  from  6q  to  6j, 
6 


-I  ‘o='P''>0.| 


2  2 

(0.  0055  rad/sec  )  t  =  0.  5  rad 


(27) 


t  =  9.  5  sec. 


(28) 


At  position  6^,  ^  will  be  maximum  so  that  D  will  begin  at  a  maxi¬ 

mum  of  about  4  ft  lb. 


When  the  viscous  damping  begins,  the  initial  conditions  are 

6  j  =  0.  7  rad 
p6j  =  -0.  1  rad/sec. 


(29) 

(30) 

The  differential  equation  which  describes  viscous  damping  phase  is 
given  by  equation  16  and  is  numerically  equal  to 

2 

p  6  +  0.  25  p6  +  0.  0156  6  =  0. 


Viscous  damping  is  defined  as  damping  which  is  porportional  to  the  rate  of  change. 


(31) 
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This  differential  equation,  with  the  initial  conditions  specified  by 
equations  29  and  30,  can  easily  be  solved  by  transforming  into  Laplacian 
space. 

2-  _ 

p  6  +  0.  25  p6  +  0.  0156  6  =  6j  p  +  (p6)j  +  0.  256^  =  0.  7p  -  0.  1  +  0.  175.(32) 

It  can  be  seen  that  the  right  side  of  Laplacian  equation  32  expresses 
the  initial  conditions.  Solving  for  6  gives 

6  =  (0.  7p  +  0.  075)/(p  +  0.  125)^  (33) 

Transforming  back  into  the  time  domain  gives 

6  =  [0.  7  -  0.  0125  (t  -9.  5)]  exp  -  [0.  125  (t  -9.  5)].  (34) 

6  is  shown  in  Fig.  8  as  a  function  of  time.  From  Fig.  8,  it  can  be 
seen  that  6  reaches  point  2  (6^  =  0.  087  rad  =  5“)  23.  5  sec  after  thrusting 
begins. 

As  b  becomes  small,  so  does  p6,  and  so  do  T  and  D.  Since  the 
valve  opening  times  will  probably  restrict  T  and  D,  consider  the  case 
where  the  minimum  T  and/or  D  is  0.  2  ft  lb.  By  the  time  this  happens, 

6^  =  0.  08  radian  (or  4.  5  degrees),  which  is  within  the  6-degree  accuracy 
required. 

At  this  point  (point  3),  it  can  be  seen  from  Fig.  8  that  pb^  is 
-0.  015  rad/sec.  At  this  point,  T  must  cut  off  so  that  only  the  damping 
torque  remains.  So  in  the  increment  3  to  4,  equation  31  reduces  to 

2 

p  6  +  0.  25  p6  +  0  =  0.  (36) 

Thus  6^  =  0.  08  radian  and  p6^  =  -0.  015  rad/ sec  are  the  initial  conditions 

of  the  increment.  The  boundary  condition  is  that  =  0.  2  ft  lb.  D  =  0.  2  ft  lb 

2  ,  2 

gives  an  acceleration  at  6^  of  p  =  -0*  00125  rad/sec  ,  which  is  equiva¬ 
lent  to  an  angular  velocity  of  p6^  =  -0.  005  rad/ sec. 

The  solution  may  be  found  by  the  Laplace  transformation 
2-  - 

p  6  +  0.  25  p6  =  initial  conditions  =  6^p  +  p6^  +  6^  (36) 

to  be 

6^  ^  =  0.  02  +  0.  06  exp  (24.  6  -  t)/4  (37) 

in  the  increment  3  to  4. 
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6^  is  that  point  at  which  D  =  0.  2  ft  lb.  This  is  found  to  occur  at 

t  =  29  sec.  This  can  be  found  by  differentiating  equation  37  and  using 
the  boundary  condition  p6^  =  -0.  005  rad/sec,  . 

p6^  =  -0.  005  =  -0.  015  exp  [(24.  6  -t)/4]  (38) 

t  =  (41n  3)  +  24.  6  =  29  sec  (39) 

so  that 

6^  =  0.  02  +  0.  06  exp  (24.  6  -t)/4  =  0.  04  rad.  (40) 

By  the  time  D  and  T  reach  their  minimum  values,  the  attitude  of  the 
vehicle  is  well  within  the  desired  limits.  From  this  point  on,  T  and  D 
can  be  programed  in  a  simple  manner  to  keep  within  the  5-degree  limits. 
For  example,  T  and  D  can  be  given  constant  values  and  applied  in  a 
manner  to  add,  subtract,  or  be  zero  by  a  simple  decision  circuit  based 
on  whether  6  times  p6  is  negative  or  positive.  The  proof  of  this  method 
was  shown  by  Dr.  R.  M.  Howe’^ 


SECTION  6  -  CONCLUSIONS 

In  order  to  control  the  attitude  of  a  meteorological  rocket  in  the  man¬ 
ner  proposed,  the  longitudinal  angular  momentum  should  be  great  enough 
for  stability,  but  small  as  compared  to  the  lateral  angular  momentum 
imparted  by  the  control  torque.  The  torque  could  be  applied  by  a  cold 
propellant  discharged  by  a  bang-bang  servo  mechanism.  The  maximum 
thrust  of  the  jet  that  would  be  required  might  be  about  2.  0  lb  if  the  jet 
were  placed  65  in.  in  front  of  the  center  of  mass. 

Based  on  the  results  of  this  study,  it  appears  that  the  design  and 
manufacture  of  a  low  cost  attitude  control  system  for  a  meteorological 
rocket  is  practical.* 


*  Development  and  Test  of  a  Simple  Attitude  Control  System  for  Small  Rockets,  Department  of 
Aeronautical  and  Astronautical  Engineering,  University  of  Michigan  Research  Institute, 

Ann  Arbor,  Michigan,  November'  1960. 
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